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ABSTRACT 

We report the results of an extragalactic point source search using the 61 
and 94 GHz (V- and W-band) temperature maps from the Wilkinson Microwave 
Anisotropy Probe (WMAP). Applying a method that cancels the "noise" due 
to the CMB anisotropy signal, we find in the |6| > 10° region 31 sources in the 
first-year maps and 64 sources in the three-year co-added maps, at a 5a level. 
The la position uncertainties are 1.6' and 1.4' each. The increased detections 
and improved positional accuracy are expected from the higher signal-to-noise 
ratio of WMAP three-year data. All sources detected in the first-year maps 
are repeatedly detected in the three-year maps, which is a strong proof of the 
consistency and reliability of this method. Among all the detections, 21 are 
new, i.e. not in the WMAP three-year point source catalog. We associate all 
but two of them with known objects. The two unidentified sources are likely to 
be variable or extended as observations through VLA0, CARMaE and ATCA§ 
all show non-detection at the nominal locations. We derive the source count 
distribution at WMAP V-band by combining our verified detections with sources 
from the WMAP three-year catalog. Assuming the effect of source clustering is 
negligible, the contribution to the power spectrum from faint sources below 0.75 
Jy is estimated to be (2.4 ±0.8) x 10~ 3 fiK 2 sr for V-band, which implies a source 
correction amplitude A = 0.012 ± 0.004 /iK 2 sr. 

Subject headings: cosmic microwave background — cosmology: observations - 
methods: data analysis 



1 The Very Large Array. See |http : //www . via . nrao"! edu/ 

2 The Combined Array for Research in Millimeter-wave Astronomy. See http://www.mmarray.org/. 
3 The Australia Telescope Compact Array. See |http : //www . narrabri . atnf . csiro . au7 
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Introduction 



The Wilkinson Microwave Anisot ropy Probe was laun ched in 2001 for making precise 
measurements of the CMB anisotropy (IBennett et al.ll2003af ). Since on small scales the CMB 
signal in the WMAP full-sky maps is primarily contaminated by microwave emission from 
extragalactic point sources, it is important to detect and mask out point-source-contaminated 
pixelgj. Accordingly, two extragalactic point source searches were mad e in WMAP first-year 
maps and three-yea r co-added mapg§ by the WMAP science team ( Bennett et al. 2003c . 



1/2 

Hinshaw et al.l 120071 ). The temperature maps were first weighed by N Q l s and then filtered 



by bi I '(bfCf ' M B -\-Cf mse ) in harmonic space, where N j, s is the number of observations, bi is the 
WMAP beam transfer function, Cf MB is the CMB angular power spectrum and C™ mse is the 
noise power. Peaks that are greater than 5a in any band of the filtered maps were interpreted 
as source detections. This procedure yielded a catalog of 208 point sources in the first-year 
data and an enlarged catalog of 323 point sources in the three-year data, with a position 
uncertainty of 4' for both searches. A point source mask of nearly 700 objects was then 
cons tructed to inc l ude a ll the 323 WMAP directly detected sourc es, along with the source s 
from lStickel et al.l (119941 ) . s ources with 22 GHz flux es > 0.5 Jy from lHirabayashi et al.l (120001 ) . 



flat-spect rum objects fromlTerasranta et al 



(200 



survey of iPerlman et al.l (119981 ) and lLandt et al.l (120011 ). Each source was masked to a radius 



), an d sources from the X-ray/radio blazar 



of 0.6°. On top of that, the power spectra of the unresolved sources were modeled into 



a 



(M0 



A (^Y C^f J l ' l>) 



with A = 0.014±0.003 /iK 2 sr, (5 = -2.0, u Q = 40.7 GHz. Here wj is the window function that 
gives the combined smoothing effec ts from the beam and the finite sky map pixel size, and 
denotes a pair of DA indices (IHinshaw et al.ll2003l . 120071 ) . This model was subtracted 
from the WMAP derived power spectra preceding any cosmology analysis. 

In this paper, we present a new extragalactic point source search in the WMAP V- and 
W-band data, using a method that has no CMB dependence. The main purpose of this work 
is to find more point sources in these two bands by lowering the noise level as only 76.2% 
and 37.5% of the WMAP three-year cataloged sources are identified with > 2a confidence 
in V- and W-band, respectively. Secondly, we are interested in seeing sources that might 



4 The Hierarchical Equal Area isoLatitude Pixelization (HEALPix) of the sphere is used to define WMAP 
map pixels on the sky in Galactic coordinates. A resolution level r =9 is chosen, which corresponds to 

2003 bi 



3,145,728 pixels with a pixel resolution of 0.115° (Be nnett et al 

5 Data products are available on-line through 



the 



LAMBDA 



website: 



http : / /lambda . gsf c . nasa . gov/ product /map. 
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peak at these high frequency bands, for instance, the so-called Gigahertz-Peaked Spectrum 
(GPS) sources which are be lieved to be at the very early state of the evolution of powerful 
radio sources (jO'Deal Il998l ). A third motiv a tion arises because in a recent re-analysis of 
WMAP temperature maps by lEriksen et al.l (120071 ). an unexpected systematic discrepancy 
was found between the V- and W-band power spectra on small ang ular scales (few percent at 
I > 300). This was partly explained by iHuffenberger et al.l (120061 ) to be an over-subtraction 
of unresolved point sources, and the significance of this difference was reduced from 3cr to 
2cr by applying a smaller source correction A = 0.011 ± 0.001 /iK 2 sr. This amplitude was 
later revis ed to A = 0.012 ± 005 u K 2 sr for t < 500 and A = 0.015 ± 0.005 /iK 2 sr for 
I > 500 in IHuffenberger et al.l (120071 ). Therefore, we also aim at getting a better constraint 
on the unresolved source contamination by using direct measurements. 



2. Methodology 

We start from the V- and W-band WMAP full sky temperature maps (the CMB dipole 
has been removed). Since WMAP o bserves the sky convolved with different beam patterns 
at different bands ( Page et al.l I2OO3I ) . we first smooth the high resolution W-band map to 



match V-band resolution. A cutoff radius 9r c of 1° is chosen for the solid angle integrations. 
A smoothing function S(9) is defined as a polynomial of an infinitely smooth function with 
compact support, 

S(9) = a nSn (9), (2) 



where 




nd ■ < 9i 



s n {0) = { 7 " ^~ T1) Rc (3) 

9 > 9 Rc . 

The coefficients a n are found by minimizing the difference between the Legendre transform 
of the smoothing function 

Si = 2tt j S{9)P l {cos9)d{cos9) (4) 

and the ratios of WMAP V- and W-band beam transfer functions 

Si = bf/bf. (5) 

We start with the 1 st power of the polynomial and truncate it after the 10 th power, so only 
the coefficients a\ to aio are non-zero, as this is sufficient to give a good fit of bf /b™ ' . The 
profile of S(9) is given in Figure [TJ The W-band map is degraded to V-band resolution by 
convolving with this smoothing function. 
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We then subtract the smoothed W-band map from the V-band map. Since the CMB 
anisotropy is from gravitational perturbations and Thomson scattering of blackbody radia- 
tion, all anisotropy-induced "noise" is removed at this step. One assumption we made here 
is that the point sources we are looking for do not have a blackbody spectrum. This is a 
reasonable assumption as the sources o bserved at WMAP fr equencies are found to be pre- 



dominately flat spectrum, i.e. F v oc u° ( jBennett et al.l 12003d ) . The V-band map, smoothed 



W-band map and the residual map are shown in Figure [2j Regions around the ecliptic poles 
are seen to be less noisy compared to other parts in the sky since these areas were heavily 
observed by WMAP. And because the Galactic emission in V-band is free-free dominated 
and in W-band is dust dominated, the residual map also reflects the ratio of free-free emis- 
sion to dust emission, i.e., red represents strong free- free emission areas and blue represents 
strong dust emission areas. 

If there is a point source with flux / in the residual map and assuming the contribution 
from the overlap of other point sources is negligible, the resulting intensity at a close- by 
pixel i can be fit into 

Vi - wf M = fbj + a. (6) 

Here Vi and Wf are pixel temperatures in the V-band and smoothed W-band maps; b\ is 
short for b v (9i), where b v (9) is the V-band beam profile normalized to 1 at 9 = and 9i is 
the angular distance from pixel i to the considered source; q is the local baseline. The x 2 °f 
the above fit is 

ij 

with 

e l = (V, - Wf M ) - (fbl + Ci) (8) 

and 

fyj = dijOy + S ik S jk (Jw- (9) 

k 

Here ay is the noise at pixel i in the V-band map and aw is the noise at pixel k in the pre- 
smoothed W-band map; S is the previously-defined smoothing function. The summation over 
k is a sum over all the pixels within 9r c of pixel i and j as the W-band noise is correlated 
after the smoothing process. A least \ 2 fit then gives 



XX 1 ( ? ) ( 1 ) ( f c ) =EVW - <") ( i ) • 



(10) 
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Defining 



M= (M 00 M 01 \ =y (b lbj bi \ 

\M W M U J \ bj 1 J' [ } 
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the intensity of the source can be written as a weighted sum of the intensities in the sur- 
rounding pixels, 



{V 3 - w. 



SM\ 



= Y^mw-w™). (12) 

i 

Computing K(9) to a cutoff radius of 1° at 32 random points in the sky, we find that the 
profiles of K(6) all share the same shape, as shown in Figure [31 Hence we adopt the same 
polynomial form for the smoothing function S(8) to fit K(6), using the collective data from 
all the 32 sky points. After filtering the residual map with K(6), we minimize the noise of 
the map and as well obtain a better presentation of the point source flux distribution of the 
V - W sky 

We then look for peaks greater than 5a in the |fe| > 10° region of the filtered map. Note 
that the WMAP scan pattern gives non-uniform observation numbers among pixels: greatest 
at the ecliptic poles, high at latitude ±45°, and least in the ecliptic plane. We thus divide 
the cut sky into 200 rings with different ecliptic latitudes and compute a separately within 
each ring. When more than one > 5a detection lay within a 3-pixel by 3-pixel area, the 
brightest one is chosen as a source detection. The position coordinates of a source are then 
carefully determined by fitting a point source profile in its 9-pixel neighborhood and locating 
the local maximum. This process results in 31 point sources in the first-year maps and 64 
point sources in the three-year co-added maps. Since all the sources found in the first-year 
maps are also detected in the three- year maps, to avoid redundancy only the three-year point 
source catalog is given (Tabled]) and will hereafter be discussed in this paper. Figure [5] offers 
an overview of the source distribution on the sky. Some galactic sources are picked in our 
search due to the loose |fe| > 10° criterion. 



3. Identification 



We cross-correlate our source list with the latest WMAP point source catalog (IHinshaw et al 
20071 ). If a source is within one half of the V-band beam width (0ff WHM ~ 10') to a cataloged 



WMAP source, we tag it. We find 43 of our sources are present in the WMAP catalog, which 
leaves us 21 new detections. 



The V-band flux densities of all our sources except for source 0540-0242 and 0734-4109 
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are calculated through 

F u = J ' I v lo{6) cos 9dQ = ^\ To J T lu(6) cos 6dtt, (13) 



wher e T = 2.725±0.002 K is the CMB monopole temperature from COBE (IMather et al. 



19991 ); B y is the Planck function; T is the V-band temperature measurement (CMB dipole 
deducted); The solid angle is integrated to a radius of 1 degree. A weighting function u(9) is 
introduced in the above equation, which is constructed with the V-band beam profile in the 
center and a negative Gaussian outer ring (shown in Figure [5]). The generated weights help 
to compensate the flux dilution caused by the big beam of WMAP and also to enhance the 
contrast of the point source flux to the background. We enforce 

u{0)<m = o, (14) 



so that any flux that spreads across the whole integration field will be ignored. For source 
0540-0242, since the bright source Orion B falls into the negative ring, using the standard 
weighting function would result in a negative flux. The same concern holds for source 0734- 
4109, as it is located in the Gum Nebula complex and is surrounded by many radio emitting 
blobs. Therefore, the V-band flux densities of these two sources are calculated using equation 
(13) wit hout applying the w eights. We use the flux densities of the 43 WMAP cataloged 



sources (IHinshaw et al.l 120071 ) as a check of our flux estimates. As seen in Figure El these 
two sets of fluxes agree very well, with our values slightly off by a factor of 1.02±0.03 in 
the mean. This is likely a consequence of the negative weights in our flux estimator. We 
thus correct our flux estimates by scaling them up. The uncertainty on the flux density 
is given by the statistical error of the fit and the rms of a simulated foreground-free CMB 
V-band map (see §4) under the same flux estimator. We then search in the NASA/IPAC 
Extragalactic Database (NED§) for radio sources within 9^j WHM to our sources and inspect 
their SEDs. If the V-band flux of our source can be reasonably fitted into the SED of a 
known radio source, we tag it. This procedure helps identify 15 more sources in our catalog 
and most of them are found to have a flat spectrum. 

For source 0358+3642 and 0402+3614, several radio sources are found around our po- 
sitions, but no one alone can explain the flux of either of our sourc es. To have a direct 



view of this source-populated region, we obtain a 2° x 2° GB6 map ([Gregory et al.lll996l ) 
centered on the relatively brighter source 0402+3614, which immediately reveals that both 
sources lay within the boundary of the diffuse nebula NGC 1499 (Figure [7J). We over-plot 



See http : //nedwww . ipac . caltech. edu/ 
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radio sources from 5 GHz surveys MG2 flLaugston et alJll990h . MG3 flGriffith et alJll990h 



and GB6 on this radio image and find that although our sources are neither aligned with 
the strong radio emission bands nor coincident with any of these 5 GHz radio sources, they 
are largely within the WMAP V-band beam size. We hence consider that our sources are 
caused by the beam smearing of WMAP and detected at the position where the composite 
flux reaches its maximum. A double check in the corresponding WMAP V-W sky map (also 
shown in Figure [7]) confirms the general shape of the nebula and verifies that our sources 
are indeed the peaks in this blended emission area. For source 1620-2536, although radio 
source PMN J1620-2529 is only 7.6' away, we cannot securely associate them as the only 
flux measurement of PMN J1620-2529 was made in 1994 (318 mJy at 4.85 GHz), and was 
one order of magnitude below the V-band flux of our source (2.2 Jy). If they were the 
same source, a power spectrum index of ~0.8 is required which indicates a negative value 
in the V-W map. Moreover, this source failed to be detected in the latest Combined Ra- 



dio All-sky Targeted Eight GHz Survey (CRATES, iHealey et al.l 120071 ). However, we know 



from its coordi nates that our sourc e resides in the p Ophiuchus nebula and a 1° x 1° PMN 



map (Figure [8j iGriffith et al.l Il994l ) further shows that it coincides with a strong emission 
region within the nebula. We double check the WMAP V-W sky map of the same field (also 
given in Figure EJ) and find the brightness distribution is in good accordance with that in 
the PMN map. Similarly, we find source 0734-4109 embedded in the Gum nebula complex, 
but no strong radio point source was reported in its immediate surroundings. Therefore, no 
one-to-one connection can be made either. 

Following the above steps, we solidly verify 58 sources in our catalog, and relate 4 
more sources to known objects without attaching them to any single point source. For 
the two source candidates which we fail to associate with any known objects, snapshots 
were taken on source 0347+0733 at VLA 3.6 cm and CARMA 3 mm bands, and on source 
1144-3902 at ATCA 3.5 cm and 6.25 cm bands during late 2006 to early 2007. None of 
the above observations show any evidence of strong radio sources existing at the nominal 
positions. This non-detection could be due to a large extent for these sources so that the 
small beams of the radio interferometers failed to catch them. Or it could be due to the 
variability of these sources. As we go through individual year WMAP data, we find that 
the V-band flux densities of these two sources did vary from year to year, given in Table 
|2j To judge the significance of this variability, we bin the map pixels (excluding those with 
|6| ^ 10°) by different observation numbers, and fit the linear relation between median 
variance in each bin and the observation numbers (see Figure [9} to obtain the flux variance 
per observation crQ median - Since the three year flux variance can be approximated by a 
X 2 with 2 degrees of freedom, the mean variance per observation can then be calculated 
as a 2 = (oq) /N b s , where (ctq) = a 2 , median / In 2. We thus find that the flux variance of 
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source 0347+0733 is 4.2 times as much as the expected level, and the variance for source 
1144-3902 is 1.1 times the expected level. In addition, as a 5cr detection threshold can cause 
a few false detections in the 3 million pixels of WMAP, these two unidentified sources could 
also be spurious. 

We find seven of our sources lying outside of the WMAP point source mask. The 5 GHz 
count erparts of our sourc es are id entified by cross-co r relating our source catalog wit h the 
GB6 dGreeory et al.l ll996h. PMN jGriffith et alill994j . Il995l . IWright et al]ll994l . Il996h . and 
Kuehr et al.l (119811 ) catalogs. We also check our source list wit h the results from two recent 
large sky radio survey: C RATES (8.4 GHz.lHealey et al.ll2007l ) and the Australia Telescope 
20 GHz Survey fAT20G. iMassardi et al.l l2007h . If a source is cross-listed in the CRATES 
catalog or the Bright Source Sample (BSS) of AT20G survey, we mark it (see Tabled]). 



Discussion 



As specified in §2, we remove all the CMB signal from our search by subtracting the 
smoothed W-band map from the V-band map. Hence, noise in the residual map comes 
primarily from the detector noise, which can be reduced as l/y/i by integrating longer. 
This trend is clearly demonstrated in Figure [inj where we plot the noise dependence on 
the observation numbers in the filtered V-W map. And we have also proved it by finding 
over \/3 times as many sources in the three-year data as in the first-year data. This is a 
distinct advantage of our method compared to the one used by the WMAP science team. 
Recall that a main step in their source finding procedure is to filter the individual band 
maps with bij '(bfCf MB + C nmse ), where the CMB signal is a "noise" that does not inte- 
grate down by repeated observations. To estimate the influence of this CMB "noise", we 
generated monopole- and dipole-free CMB maps for each WMAP band using cosmological 
parameters f avored by the c ombined data sets of WMAP 3 year data, 2dF Galaxy Red- 



shift Survey (I Cole et a 



vey 



2005), BOOMERanG and ACBAR fMontrov et al.ll2006l. iKuo et al. 



20041), CBI and VSA flReadhead et all 12004 1. Dickinson et al 



(ITegmark et al.l 120041 . lEisenstein et al. 



20041 ) . Sloan Digi tal Sky Sur- 



2005h. Super nova Legacy Survey ( Astier et al. 



20061 ) and Supernova "Gold Sample" (IRiess et al.l 120041 ) . assu ming a ACDM model . We 
then repeat the same procedure of weighting and filtering that iHinshaw et al.l (120071 ) did, 
and find the rms fluctuation level in the resultant maps to be 0.27 Jy, 0.41 Jy, 0.36 Jy, 0.27 
Jy and 0.14 Jy for K-band through W-band, individually. Considering that the detected 
sources are mostly about 1 Jy, this intrinsic CMB "noise" may have caused some missing 
sources at a 5cr level. 



In addition, since the V- and W-band maps have comparatively higher angular resolution 
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among the five WMAP bands, we achieve a higher positional accuracy in our search. The 
la position uncertainty is 1.6' and 1.4' for the first-year and three-year result, respectively. 
We also achieve a low detection threshold of 0.7 Jy at V-band in our search, attributed to 
the elimination of the CMB signal and the relatively low foreground contamination in the 
V- and W-band sky. However, it is noticed that we did not detect all the sources with 
Fy greater than 0.7 Jy in the WMAP three- year point source catalog. This is because the 
above-mentioned subtraction step not only removed the CMB signal but also reduced the 
intensity of possible sources. Assuming a flat spectrum, the W-band temperature of a source 
is about 1/2 the V-band temperature. Therefore some strong sources in V-band could be 
diminished in the residual map. 



Besides the searches done by the WMAP science team, iNie fc Zhang! (120071 ) andlLopez-Caniego et al. 



( 2007 ) independently carried out a point source search in the WMAP map s. INie &: Zhang 
(120071 ) applied the cross-correlation method developed by iHermsenl (119831 ) on the WMAP 
first year Q-band m ap and found 20 new s ource s at |6| > 10°, out of which we found 9. In the 
non-blind search by lLopez-Caniego et al.l (120071 ). they first generated an input catalog of 2491 



sourc es from the 4.85 GHz GB6 and P MN surveys, 1.4 G Hz NVSS survey (jCondon et al. 



1998|) and 0.8 43 GHz SUMSS survey (IMauch et all 12003). and then applied the MHW2 
wavelet filter ( Gonzalez-Nuevo et al. 20061 ) to 14.6° x 14.6° patches surrounding the sources 
and picked the 5a detections. 15 new sources were found at V-band with |6| > 10°, out of 
which we also found 9. These cross-detections, along with the 43 WMAP cataloged sources, 
prove that our approach is both competitive with and complementary to the other currently 
adopted methods for detecting point sources in the WMAP data. Furthermore, our V-band 
flux estimates for these common sources are in general consistent with the estimates from 
other methods, which verifies the effectiveness of our flux estimator. 

Combining our identified detections (excluding galactic sources) and the cataloged 
WMAP sources, we modeled the V-band source count distribution dN/dS as a power law 
k(S/[1 Jy])~^, with k = 17.6 ± 3.3 and (3 = 2.4 ± 0.1. For Poisson-distributed sources, the 
contribution to the power spectrum by unresolved sources can be estimated by 



C 



9\y) 



dS——S 2 /iK 2 sr, 
db 



(15) 



where S c is the flux limit and g(y) = (dB^/dT)^ 1 . Integrating this model with S c = 0.75 Jy 
resulted in a contamination level of Cy = (2.4 ± 0.8) x 10~ 3 /iK 2 sr for V-band. Fitting this 
value to the power spectrum model of the unresolved sources, equation (1), we obtained a 
Q-ba nd amplitude o f A = 0.012 ± 0.004 nK 2 sr, wh i ch is consistent with the estimates from 



both iHinshaw et all (120071 ) and iHuffenberger et al.l (120071 ) . 
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Conclusions 



We present a CMB independent method for point source detection and apply it to the 
WMAP first-year and three-year maps. The main result of this search is a point source cata- 
log of 64 members, out of which 21 are newly detected in WMAP data. Two sources remain 
unidentified. We are not sure whether they are variable or extended sources. More years 
WMAP data are needed to further identify them. We demonstrate that our method is both 
competitive with and complementary to other currently adopted methods of point source 
searching in the WMAP data, with a low detection threshold as well as a good positional 
accuracy (~ 15% of the W-band FWHM). And our source list can be a good supplement to 
the existing WMAP point source catalog and the point source mask. Contamination from 
unresolved point sources is at a negligible level of (2.4 ±0.8) x 10 -3 fiK 2 sr at WMAP V-band, 
which implies a correction value A = 0.012 ± 0.004 fiK 2 sr. 

For most sources, with F v oc u°, the K, Ka and Q bands have a higher signal to noise 
ratio and therefore c a n yield more s ource dete ctions as proven in the sea rch results from 
Hinshaw et al.1 (120071 ). iNie fc Zhang! (120071 ) and lLopez-Caniego et al.l (120071 ). However, part 
of the intention of this work was to look for interesting sources that might peak around the 
V and W bands. Nevertheless, among all the identified sources in our catalog, we find 39 
quasars, 10 galaxies, 9 HII regions and 4 other radio sources. There is no sign of a significant 
different class of new objects. 

We are planning to extend this work by including the WMAP Q-band data. We an- 
ticipate having more source detections so that a more accurately model of the source count 
distribution and furthermore a better estimate of the contribution from unresolved point 
sources to the CMB power spectrum can be reached. 
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Table 1. Point Source Catalog -Three Years 



RA 
hms 


DEC 
dms 


WMAP ID 


Type 


F V 

[Jy] 


5GHz ID 


Identified 
/Masked 3 


Notes b 





00 06 


14 


-06 23 44 


WMAP 


J0006-0622 


G 


2.2+0.3 


PMN J0006-0623 


Y 


/ 


Y 


t 


01 08 


47 


01 37 30 


WMAP 


J0108+0135 


QSO 


2.4±0.3 


GB6 J0108+0135 


Y 


/ 


Y 


t 


01 36 


58 


47 51 27 


WMAP 


J0137+4752 


QSO 


3.5±0.3 


GB6 J0136+4751 


Y 


/ 


Y 


i 


02 10 


46 


-51 00 07 


WMAP 


J0210-5100 


QSO 


2.7±0.3 


PMN J0210-5101 


Y 


/ 


Y 


t, O 


02 37 


51 


28 45 47 


WMAP 


J0237+2848 


QSO 


2.7±0.3 


GB6 J0237+2848 


Y 


/ 


Y 


t 


03 19 


52 


41 30 35 


WMAP 


J0319+4131 


G 


4.9±0.4 


GB6 J0319+4130 


Y 


/ 


Y 


t, Per A 


03 21 


53 


-37 05 50 






RadioS 


1.9±0.3 


PMN J0321-3711 c 


Y 


/ 


Y 


*, For A west lobe 


03 47 


16 


07 34 00 








0.7±0.3 




N 


/ 


N 




03 58 


38 


36 41 53 






HII 


2.4±0.3 




Y 


/ 


Y 


In NGC 1499 


04 02 


41 


36 15 06 






HII 


4.6±0.4 




Y 


/ 


Y 


In NGC 1499 


04 03 


55 


-36 05 25 


WMAP 


J0403-3604 


QSO 


4.6±0.4 


PMN J0403-3605 


Y 


/ 


Y 


t, o 


04 23 


20 


-01 19 16 


WMAP 


J0423-0120 


QSO 


9.2±0.5 


PMN J0423-0120 


Y 


/ 


Y 


t 


04 40 


39 


-43 36 20 


WMAP 


J0440-4333 


QSO 


2.4±0.3 


PMN J0440-4332 


Y 


/ 


Y 


t, o 


04 55 


45 


-46 17 22 


WMAP 


J0455-4617 


QSO 


3.6±0.3 


PMN J0455-4616 


Y 


/ 


Y 


t, o 


04 57 


15 


-23 23 34 


WMAP 


J0457-2322 


QSO 


1.9±0.3 


PMN J0457-2324 


Y 


/ 


Y 


t, o 


04 57 


16 


-66 25 31 






HII 


3.1±0.3 


PMN J0456-6624 


Y 


/ 


Y 


f, in LMC 


05 17 


10 


-69 20 16 






HII 


1.0±0.3 


PMN J0518-6914 


Y 


/ 


N 


f, in LMC 


05 19 


55 


-45 48 59 


WMAP 


J0519-4546 


G 


3.1±0.3 


Uy 0518-45 


Y 


/ 


Y 


o, Pic A 


05 22 


18 


-68 01 44 






HII 


1.7±0.3 


PMN J0522-6757 


Y 


/ 


N 


f, t, in LMC 


05 23 


22 


-36 26 25 


WMAP 


J0523-3627 


G 


3.0+0.3 


PMN J0522-3628 


Y 


/ 


Y 


o 


05 35 


16 


-05 23 29 






HII 


290.8±8.0 




Y 


/ 


Y 


Ori A 


05 38 


25 


-69 07 37 






HII 


26.2±0.9 


PMN J0538-6905 


Y 


/ 


Y 


f, f, in LMC 


05 38 


35 


-44 06 17 


WMAP 


J0538-4405 


QSO 


5.2±0.4 


PMN J0538-4405 


Y 


/ 


Y 


t: O 


05 40 


33 


-02 42 15 






RadioS 


8.1±0.4 


PMN J0541-0239 c 


Y 


/ 


N 


t 


05 41 


43 


-01 53 40 






HII 


45.6±1.4 


PMN J0541-0154 


Y 


/ 


Y 


t 


06 07 


58 


-06 25 52 






HII 


6.9±0.4 


PMN J0607-0623 


Y 


/ 


Y 


t, * 


06 09 


42 


-15 40 58 


WMAP 


J0609-1541 


QSO 


2.2±0.3 


PMN J0609-1542 


Y 


/ 


Y 


t, o 


06 35 


17 


-75 16 21 


WMAP 


J0635-7517 


QSO 


3.0±0.3 


PMN J0635-7516 


Y 


/ 


Y 


t, o 


07 34 


13 


-41 09 31 






RadioS 


3.9±0.3 




Y 


/ 


N 


In Gum Nebula 


08 36 


25 


-20 16 52 


WMAP 


J0836-2015 


QSO 


2.1±0.3 


PMN J0836-2017 


Y 


/ 


Y 


t, o 


08 41 


00 


70 54 22 


WMAP 


J0841+7053 


QSO 


1.7±0.3 


GB6 J0841+7053 


Y 


/ 


Y 


t 


08 55 


05 


20 05 26 


WMAP 


J0854+2005 


QSO 


4.3±0.4 


Uy 0851+20 


Y 


/ 


Y 


t 


09 09 


09 


01 26 01 


WMAP 


J0909+0118 


QSO 


1.4±0.3 


GB6 J0909+0121 


Y 


/ 


Y 


i 


09 27 


04 


39 04 43 


WMAP 


J0927+3901 


QSO 


4.0±0.3 


GB6 J0927+3902 


Y 


/ 


Y 


t 


10 58 


32 


01 33 17 


WMAP 


J1058+0134 


QSO 


4.0±0.3 


Uy 1055+01 


Y 


/ 


Y 


i 


10 58 


54 


-80 03 13 


WMAP 


J 1059-8003 


QSO 


2.4±0.3 


PMN J1058-8003 


Y 


/ 


Y 


t, o 


11 44 


20 


-39 02 18 








0.8±0.3 




N 


/ 


N 




12 29 


09 


02 02 59 


WMAP 


J1229+0203 


QSO 


13.1±0.6 


GB6 J1229+0202 


Y 


/ 


Y 


t 


12 30 


49 


12 22 51 


WMAP 


J1230+1223 


G 


9.4±0.5 


GB6 J1230+1223 


Y 


/ 


Y 


Vir A 


12 56 


11 


-05 47 31 


WMAP 


J1256-0547 


QSO 


16.6±0.6 


PMN J1256-0547 


Y 


/ 


Y 


i 


13 10 


49 


32 23 40 


WMAP 


J1310+3221 


QSO 


2.3±0.3 


GB6 J1310+3220 


Y 


/ 


Y 


t 


13 25 


33 


-42 59 19 






G 


25.7±0.9 


PMN J1325-4257 c 


Y 


/ 


Y 


f, *, o, Cen A 


13 37 


43 


-12 57 15 


WMAP 


J1337-1257 


QSO 


6.4±0.4 


PMN J1337-1257 


Y 


/ 


Y 


i 


16 13 


42 


34 12 38 


WMAP 


J1613+3412 


QSO 


2.7±0.3 


GB6 J1613+3412 


Y 


/ 


Y 


t 


16 20 


23 


-25 36 11 






RadioS 


2.2±0.3 




Y 


/ 


N 


In p Oph nebula 
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Table 1 — Continued 



RA 
hms 


DEC 
dms 


WMAP ID 


Type 




F V 

[Jy] 




5GHz ID 


Identified 
/Masked a 




Notes b 


16 


35 


03 


38 07 91 


WMAP 


J1635+3807 


QSO 


4 


7dz0 


4 


OR6 T1 635^3808 


Y 


/ 


Y 


+ 
+ 






16 


42 


57 


39 48 38 
oy oo 


WMAP 


J1642+3948 


QSO 


5 


4-J-Q 


4 




Y 


/ 


Y 


t 






17 


0.3 


31 


-62 10 36 

1U ou 


WMAP 


J1703-6213 


G 


2 


,1±0. 


.3 


PMN 11703-6212 


Y 


/ 


Y 


+ 


o 




17 


20 


30 


-00 59 15 






G 


2 


.6±0. 


.3 


PMN J1720-0058 


Y 


/ 


Y 


t, 


* 




17 


33 


00 


-13 06 22 






QSO 


4 


2±0. 


.4 


PMN J1733-1304 


Y 


/ 


Y 


*J 


t 




17 


13 


57 


-03 48 35 






QSO 


1 


,7±0. 


.4 


PMN J1743-0350 


Y 


/ 


Y 


*, 


i 




17 


51 


41 


09 41 10 






QSO 


3 


,8±0. 


.3 


GB6 J1751+0938 


Y 


/ 


Y 


t, 


*, 


t 


17 


53 


43 


28 47 04 


WMAP 


J 1753+2848 


QSO 


2 


3±0. 


3 


GB6 J1753+2847 


Y 


/ 


Y 


t 






18 


06 


38 


69 50 48 


WMAP 


J1806+6949 


G 


1 


,2±0. 


.3 


GB6 J1806+6949 


Y 


/ 


Y 


t 






19 


21 


50 


-29 14 38 






QSO 


9 


.9±0. 


.5 


PMN J1924-2914 


Y 


/ 


Y 


t, 


*, 


t, o 


19 


27 


39 


73 57 42 


WMAP 


J1927+7357 


QSO 


3 


,1±0. 


.3 


GB6 J1927+7357 


Y 


/ 


Y 


t 






21 


48 


10 


06 57 04 


WMAP 


J2148+0657 


QSO 


6 


,4±0. 


.4 


GB6 J2148+0657 


Y 


/ 


Y 


t 






21 


57 


11 


-69 41 18 


WMAP 


J2157-6942 


G 


2 


,2±0. 


.3 


PMN J2157-6941 


Y 


/ 


Y 


o 






22 


02 


41 


42 17 11 


WMAP 


J2202+4217 


QSO 


3 


0±0. 


.3 


GB6 J2202+4216 


Y 


/ 


Y 


* 






22 


18 


59 


-03 33 18 


WMAP 


J2218-0335 


QSO 


1 


,7±0. 


.3 


PMN J2218-0335 


Y 


/ 


Y 


i 






22 


25 


11 


-04 58 56 


WMAP 


J2225-0456 


QSO 


3 


,2±0. 


3 


PMN J2225-0457 


Y 


/ 


Y 


t 






22 


29 


56 


-08 33 29 


WMAP 


J2229-0832 


QSO 


3 


3±0. 


3 


PMN J2229-0832 


Y 


/ 


Y 


t 






22 


53 


57 


16 09 07 


WMAP 


J2254+1608 


QSO 


6 


,7±0. 


,4 


GB6 J2253+1608 


Y 


/ 


Y 


t 






22 


58 


09 


-27 59 58 


WMAP 


J2258-2757 


QSO 


6 


,2±0. 


A 


PMN J2258-2758 


Y 


/ 


Y 


t, 








a Three-year WMAP point source mask is considered here. 

b f and * indicate the new point sources cross-detected in lNie &~Z hang (20 071') and | L 6pcz-Canic go et all ll2007t ). | 
and o indicate the source is included in the CRATES catalog dHealev et alj 1200711 and the AT20G BSS cata- 
log (JlVlassardi et al. 2007), respectively. 

c Indicates the source has multiple possible 5 GHz identifications. The brightest one is given here. 



Table 2. Single Year V-band Fluxes of the Two Unidentified Sources 



RA 


DEC 




F yr2 


pyrs 


4 


(cr%)/N obs 


hms 


dms 


[Jy] 


[Jy] 


[Jy] 


[Jy 2 ] 


[Jy 2 ] 



03 47 16 07 33 25 1.9±0.5 1.2±0.5 0.5±0.6 0.50 0.12 
11 44 20 -39 02 18 1.5±0.5 1.2±0.5 0.9±0.6 0.08 0.07 
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0.000 0.005 0.010 0.015 0.020 
6 (radians) 



Fig. 1. — Smoothing function S(9), used to degrade the W-band map to V-band resolution. 
The cutoff radius 6 Rc is l°(f« 0.0 f 7 radians). 



-0.20 



0.20 



Fig. 2. — (a) The V-band sky map. (b) The W-band sky map smoothed to V-band resolution, 
(c) The residual sky map after subtracting (b) from (a). Regions around the ecliptic poles 
and latitude ± 45° are heavily observed by WMAP and thus appear less noisy. The maps 
are shown in the Mollweide projection in units of CMB thermodynamic temperatures (mK). 
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0.3 F 




-0.1 F .... i .... i .... i .... 1 
0.000 0.005 0.010 0.015 0.020 
8 (radians) 

Fig. 3. — The filter we used to minimize the noise and recover the point source flux distri- 
bution in the V-W residual map. Dots are superimposition of generated data at 32 random 
sky locations; Solid line is our fit. The scale of K(0) is determined by the normalization of 
the V-band beam profile as indicated by Equation (6) and Equation (12). 



+90 




Fig. 4. — An overview of the microwave sky. The pink and violet shaded regions indicate the KpO, Kp2 diffuse emission 
masks. The blue dots give the point source mask. The red and green dots together give the sources detected in our 
search; Green means it is included in the existed WMAP point source mask while red means not included. 
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-t— ' 

0> 




-0.5 1 . ... i .... i .... i .. . , 
0.000 0.005 0.010 0.015 0.020 
6 (radians) 

Fig. 5. — The weighting function (solid line) we introduced in our V-band flux estimator, 
which is constructed with a V-band beam matched central peak (dotted line) and a negative 
Gaussian (dashed line) outer ring. Application of this weighting function provides a flux 
contrast enhancement of the point source from the background. 




F v (Jy), this work 



Fig. 6. — Comparison of the point source V-band flux densities we initially computed in 
this work and those from the WMAP three-year release. The diamonds are the 43 WMAP 
cataloged point sources we rediscovered in this work. The linear fit has a slope of 1.02. 
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Fig. 7.— Left: 2° x 2° GB6 map including our source 0402+3614 and 0358+3642. White 
circle covers a 3a error area. Cross points give MG2 and MG3 sources within one half the 
WMAP V-band beam width. The X point is the closest GB6 source to our sources, which is 
13.3' away from source 0402+3614. Right: Corresponding 2° x 2° field in the filtered V-W 
map. 




Fig. 8. — Left: 1° x 1° PMN map centered on our source 1620-2536, marked by a 3a error 
circle. Cross point indicates the PMN source within one half the WMAP V-band beam 
width. Right: Corresponding 1° x 1° field in the filtered V-W map. 
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0.25 




Fig. 9. — The three years V-band flux variance dependence on observation numbers. The 
diamonds are the median sample variances measured from each bin of map elements. The 
solid line is a linear fit to the data. 



0.008 




Fig. 10. — Noise dependence on numbers of observation in the filtered V-W map. The dia- 
monds give the rms in each of the 200 rings at different ecliptic latitudes, and accordingly, 
different observation numbers. Since all the convolutions we did to the maps are mathe- 
matically linear combinations, the linear relation between a and N^J 2 is preserved, and the 
solid line gives the fit. 



